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ABSTRACT 

We use a semi-analytical model to study the impact of reionization, and the associated 
radiative feedback, on galaxy formation. Two feedback models have been considered: 
(i) a standard prescription, according to which star formation is totally suppressed 
in galaxies with circular velocity below a critical threshold (model CF06) and (ii) 
a characterization based on the filtering scale (model GOO), allowing for a gradual 
reduction of the gas available for star formation in low-mass galaxies. In model CF06 
reionization starts at z $ 15 — 20, is 85% complete by z ~ 10; at the same z, the ionized 

fraction is 16% in model GOO. The models match SDSS constraints on the evolution of 
the neutral hydrogen fraction at z < 7, but predict different Thomson optical depths, 
t c = 0.1017 (CF06), and 0.0631 (GOO); such values are within la of the WMAP 
3-yr determination. Both models are in remarkable good agreement with additional 
existing data (evolution of Lyman-limit systems, cosmic star formation history, high- z 
galaxy counts, IGM thermal history), which therefore cannot be used to discriminate 
among different feedback models. Deviations among radiative feedback prescriptions 
emerge when considering the expected HI 21 cm background signal, where a ~ 15 mK 
absorption feature in the range 75-100 MHz is present in model GOO and a global shift 
of the emission feature preceding reionization towards larger frequencies occurs in 
the same model. Single dish observations with existing or forthcoming low-frequency 
radio telescopes can achieve mK sensitivity, allowing the identification of these features 
provided that foregrounds can be accurately subtracted. 

Key words: Cosmology: theory - galaxies: formation - intergalactic medium - diffuse 
radiation 



1 INTRODUCTION 

The process of cosmic reionization is of primary importance 
in the evolution of the Universe. Reionization is a conse- 
quence of the formation of the first luminous sources that 
shine after the Dark Ages, but it also has a dramatic impact 
on subsequent galaxy evolution. Despite the recent progress 
in observational and theoretical studies, fundamental ques- 
tions such as what are the dominant sources of reionization 
and how long did the process take still require complete and 
definite answers. 

The latest analysis of Lya absorption in the spec- 
tra of the 19 highest redshift Sloan Digital Sky Sur- 
vey (SDSS) quasars (QSOs) shows a strong evolution of 
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the Gunn-Peterson Lya opacity at z ~ 6 <|Fan et al.l 
l2006t iGallerani, Choudhury fc Ferraral I2006T) . This result, 
together with the downward revision of the electron scat- 
tering optical depth to r e = 0.09 ± 0.03 in the release of 
the 3- yr Wilkinson Microwave Anisotropy P robe ( WMAP) 
data jPage et all 120071 : ISpergel et all l2007r i. is consistent 
with "minimal reionization models" which do not require 
the p resence of very massive (M > IOOMq) Populatio n III 
stars l|Choudhurv fc Ferraral [20061 : iGnedin fc Fanll2006l ). 

Choudhury & Ferrara (2005, 2006) have developed a 
self-consistent formalism which allows to jointly study cos- 
mic reionization and the thermal history of the intergalac- 
tic medium (IGM). Their semi-analytic model accounts for 
inhomogeneous IGM density distribution, three different 
classes of ionizing photons (Population III, Population II 
stars and QSOs), and radiative feedback inhibiting star for- 
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mation in low-mass galaxies. Note that Population III (Pop 
III) stars in this model are characterized by a standard 
Salpeter Initial Mass Function (IMF) extending in the range 
1 — 100 Mq. The model free parameters are constrained by 
a wide range of observational data. From this analysis, it 
emerges that hydrogen reionization starts around z ~ 15 
driven by the contribution of Pop III stars and it is 80% 
complete by z ft! 10; below this redshift, the contribution of 
Pop III stars decreases because of the combined action of 
radiative and chemical feedback. As a consequence, reion- 
ization is extended consider ably, completing only at z ~ 6 
ijChoudhurv fc Ferraral 120061 ). 

These results are consistent with the study of 
iGnedin fc Far] { 2006), who used a set of numerical simu- 
lations with different spatial and mass resolutions to make a 
detailed comparison with SDSS and WMAP data. Assuming 
the sources of ionizing photons to be Population II (Pop II) 
stars and QSOs, they find that the simulations can match 
the SDSS observations in the range 5 < z < 6.2, but do not 
have enough resolution to resolve the earliest stages of star 
formation. Depending on the fractional contribution of Pop 
III stars, the optical depth to Thomson scattering varies in 
the range 0.06 < r e < 0.10, consistent with the 3-yr WMAP 
results. 

Since it appears that minimal reionization models pro- 
vide the best theoretical frameworks for the existing data, 
we can use these models to explore the effects of reionization 
on galaxy formation, to which we will refer to as "radiative 
feedback" . 

It is well known that the temperature increase of 
the cosmic gas in ionized regions leads to a dramatic 
suppression of the formation of low-mass galaxies. Such 
suppression has been found to become effective in dark 
matter halos with circular velocity below a critical one, 
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depending on the assumed intensity of the UV background, 
the inclusion of radiative transfer effects, the numerical 
scheme adopted, and the investigated redshift range (see 
ICiardi fc Ferrarall2005l for an extensive review). This radia- 
tive feedback prescription is therefore unable to quantify the 
decrease in the amount of gas available for star formation in 
low-mass halos: it simply states that galaxies can form stars 
unimpeded provided that their halos have circular velocities 

> Vcrit- 

A perhaps more complete characterization of r adiative 
feedba ck on low-mass galaxy has been obtained bv iGnedlnl 
(2000). Based on cosmological simulations of reionization, 
this study shows that the effect of photoionization is con- 
trolled by a single mass scale in both the linear and non 
linear regimes. The gas fraction within dark matter halos at 
any given moment is fully specified by the current filtering 
mass, which directly corresponds to the length scale over 
which baryonic perturbations are smoothed in linear theory. 
The results of this study provide a quantitative description 
of radiative feedback, independently of whether this is phys- 
ically associated to photoevaporative flows or due to accre- 
tion suppression. 

In the present paper, we show that these two differ- 
ent feedback prescriptions have important effects on the 



reionization history and on the properties of the sources 
which dominate the process. Using the semi-analytic model 
of Choudhury & Ferrara (2005, 2006), we compare the re- 
sults obtained for these two competitive feedback models 
with a wide range of observational data, ranging from the 
redshift evolution of Lyman-limit absorption systems, the 
Gunn-Peterson and electron scattering optical depths, the 
cosmic star formation history, and number counts of high- z 
sources in the NICMOS Hubble Ultra Deep Field (HUDF). 
We find that in spite of this demanding benchmark of obser- 
vations, existing data are unable to discriminate among the 
two reionization histories. We therefore explore an alterna- 
tive method to break these degeneracies using future 21 cm 
experiments such as the Low Frequency Array (LOFAR) , the 
21-Centimeter Array (21CMA), the Mileura Wide Field Ar- 
ray (MWA), and the Square Kilometer Array (SKA); these 
instruments are expected to reach the sensitivity required to 
perform accurate maps of the neutral hydrogen distribution 
from the Dark Ages to the latest stages of the reionization 
process. 

The paper is organized as follows: in Section[5]we give a 
brief description of the semi-analytic formalism and present 
the results of the two radiative feedback models, comparing 
these with existing observations; in Section[3]we compute the 
predicted global 21 cm background in the two models and 
discuss its detectability with on-going and future facilities. 
Finally, Section U summarizes our conclusions. 

Throughout this paper, we adopt a flat ACDM cosmo- 
logical model consistent with 3-yr WMAP data (Spergel et 
al. 2007), with matter and cosmological constant density 
parameters Q m — 0.24 and Q,a = 0.76, reduced Hubble con- 
stant h = 0.73, baryon density Q,bh 2 = 0.022, density con- 
trast erg = 0.74, and adiabatic scalar perturbations (with- 
out running) with spectral index n s = 0.95. We also as- 
sume a Cos mic Microwave Bac kground (CMB) temperature 
of 2.725 K (|Mather et al.ll 19991 ). 



2 REIONIZATION MODELS 

In this Section we first summarize th e main features of the 
semi-a nalytical model developed by I Choudhury fc Ferraral 
J2005I) with the a dditio nal physics introduced in 
IChoudhurv fc Ferraral (|2006h . We then describe the two 
alternative radiative feedback prescriptions and compare 
the resulting reionization histories with existing data. 



2.1 Model description 

The main features of the formalism developed by Choudhury 
& Ferrara (2005, 2006) can be summarized as follows: 

• Inhomogeneous reionization: the model accounts for 
IGM inhomogeneities by adopting the procedure of Miralda- 
Escude, Haehnelt & Rees (2000). The overdensity distribu- 
tion is assumed to be lognormal. The distribution deter- 
mines the mean free path of photons, 



A m fp(z) 
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(1) 



[1 - FV(z)] 2 /3 

where Fv is the volume fraction of ionized regions and Ao 
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Table 1. Best-fit parameters for models CF06 and GOO shown in Figs. [2] and [3] For the same 
models, we also report the residual volume-averaged neutral hydrogen fraction at redshift 
6, xjji(6), the volume- averaged electron fraction at redshfit 10, rr o (10), and the Thomson 
scattering optical depth r e i (see text). 



Model 


e*,II 


/"esc, 1 1 




/esc, III 


xhi(6) 


ie(10) 


Tel 


CF06 


0.10 


0.01 


0.03 


0.68 


4 x 10~ 4 


0.85 


0.1017 


GOO 


0.10 


0.01 


0.01 


0.12 


3.9 x 10~ 4 


0.16 


0.0631 



is a normalization constant fixed by comparing with low- 
redshift observations of Lyman-limit systems. 

• Sources of ionizing photons: the IGM is treated as a 
multiphase medium, following the thermal and ionization 
histories of neutral, HII, and Helll regions simultaneously. 
Three sources have been assumed to contribute to the ion- 
izing flux: (i) Pop III stars, assumed to be distributed ac- 
cording to a standard Salpeter IMF, as suggested by the 
combination of co nstraints from source counts at z ~ 10 
and WMAP data (|Schneider et al.ll2006t ); (ii) Pop II stars, 
with Z = 0.2Z© and Salpeter IMF; (hi) QSOs which are 
significant sources of hard photons at z ^ 6. The emission 
is mo delled using the ste ll ar po pulation synthesis templates 
from Bruzual fc Chariot! (120031 ) for Pop II stars and from 
ISchaererl (|2002h for Pop HI stars. 

• Chemical feedback: the transition from Pop III to Pop 
II stars is controlled by chemical feedback and occurs over a 
prolonged epoch rather than at a precise transition redshift. 
Using the merger-tree model developed by Schneider et al. 
(2006), at each redshift we classify star forming halos as 
hosting Pop II (Pop III) stars depending on whether the 
halo itself or any of its progenitors have (have not) already 
experienced an episode of star formation^]. The fraction of 
Pop III halos decreases with time, being w 0.4, 0.32, and 0.23 
at 2 = 15, 10, and 5, respectively. At each redshift, Pop III 
stars are confined to form in halos with masses 10 s - 1O 9 M 
which are large enough to form stars but small enough to 
be relatively unpolluted by their progenitors. 

• Escape fractions: to reduce the number of model free 
parameters, Choudhury & Ferrara (2006) use a physical ar- 
gument to relate the escape fractions of Pop III and Pop 
II ionizing photons. These scale according to the number 
of ionizing photons produced and depend on a single free- 
parameter, TVabs, which represents the number of ionizing 
photons absorbed within a star-forming halo, through the 
relation 

Aabs / \ 

77 CSC = . (2) 

e*,niV 7 ,ii 

The escape fractions for Pop II and Pop III stars are then 
given by the following expressions, 



/cscJI = 1 - Min[l, 77esc] 
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where Ny,n (iV 7) iii) are the number of photons produced by 
Pop II (Pop III) stars per unit mass of stars formed, and e*,n 
(e»,in) are the Pop II (Pop III) star formation efficiencies. 

Further details on the models can be found in the orig- 
inal papers (Choudhury & Ferrara 2005, 2006). 

2.2 Radiative feedback 

A variety of feedback mechanisms can suppress star for- 
mation in mini-halos, i.e. halos with virial temperatures 
< 10 4 K, particularly if t heir clustering is taken into ac- 
count (|Kramer et al.|[2006l ). We therefore assume that stars 
can form in halos down to a virial temperature of 10 4 K, 
consistent with the interpretation of the 3-yr WMAP data 
(Haiman & Bryan 2006; but see also Alvarez et al. 2006). 
We then assume that these halos can be affected by radiative 
feedback according to two alternative prescriptions: 

(i) Following IChoudhurv fc Ferrara! l|2006l ). we assume 
that in photoionized regions halos can form stars only if 
their circular velocity exceeds the critical value 



Wcrit 



2k B T 



(5) 



where fi is the mean molecular weight, m p is the proton 
mass, and T is the average temperature of ionized regions, 
which we can compute self-consistently from the multiphase 
IGM model. Typically, for a temperature of 3 x 10 4 K, 
the minimum circular velocity is w 30 km/s, within the 
range of values found in the literature (see Section [1] and 
ICiardi fc Ferraral 120051 for a complete reference list). Note 
that the above value for v C nt evolves according to the gas 
temperature and it is not fixed to a particular value. Here- 
after, we shall r efer to this mo del as CF06. 

(ii) Following iGrredinl l|2000l) , we assume that the average 
baryonic mass Mb within halos in photoionized regions is a 
fraction of the universal value = fij,/fl m , given by the 
following fitting formula, 



Mb 
M 



h 



[l + (2 1 /3 _ l)M c /M]3 : 



(6) 



where M is the total halo mass, and Mc is the total mass of 
halos that on average retain 50% of their gas mass. A good 
approximation for the characteristic mass Mc is given by 
the linear-theory filtering mass, 



1 Since Pop III stars are assumed to be distributed according to 
a Salpeter IMF, star formation in Pop III halos always leads to 
metal-enrichment as a consequence of type-II supernova explo- 
sions. This is equivalent to the strong feedback case, fsn — 1} in 
the formalism of Schneider et al. (2006). 
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where a is the cosmic scale factor, 
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Figure 1. Rcdshift evolution of the filtering mass for the GOO 
model (dashed line), the minimum mass of star forming halos 
(dotted line), and the critical mass in ionized regions for the CF06 
model (solid line). 



47T . 



Gp 



3/2 



is the Jeans mass, p is the average total mass density of 
the Universe, and c s is the gas sound speed. Since the gas 
density and temperature are not uniform, the sound speed 
is computed according to 

2 5fc s (T) v 

c. — ■ — ■ — 

3 pirip 

where (T)v in dicates the volume- averaged gas temperature 
l|Gnedinll2000l ). Hereafter, we will refer to this model as GOO. 

In Fig. [1] the filtering mass for the GOO model (dashed 
line) and the critical mass in ionized regions for the CF06 
model (solid line) are plotted as a function of redshift. While 
these masses play different roles in the two feedback models, 
it is interesting to compare them with the minimum halo 
mass allowed to form stars (dotted line). 

2.3 Results 

The semi-analytical model depends only on four free pa- 
rameters: the star formation efficiencies of Pop II and Pop 
III stars, the parameter J7 CSC which is related to the escape 
fraction of ionizing photons emitted by Pop II and Pop III 
stars (see equation [2} , and the normalization of the pho- 
ton mean free path, Ao (see equation [1]) , which is fixed to 
reproduce low-redshift obser vations of Lyman-limit systems 
l|Choudhurv fc Ferrarall2005l ). In Figs. [2] and [3] we show how 
the two models compare to existing observational data. The 
best-fitting parameters for each model are reported in Ta- 
ble 1. Both models show a remarkable agreement with a vari- 
ety of observations. However, the two feedback prescriptions 



have a noticeable impact on the overall reionization history 
and the relative contribution of different ionizing sources. 
The main reason is that, although the two models predict 
similar global star formation histories dominated by Pop 
II stars (panels b), the Pop III star formation rates have 
markedly different redshift evolution. In fact, chemical feed- 
back forces Pop III stars to live preferentially in the smallest, 
quasi-unpolluted halos (virial temperature ^ 10 4 K) , which 

are those most affected by radiative feedback (see Figure [1]). 

In model CF06, where star formation is totally sup- 
pressed below fcrit, Pop III stars disappear at z ~ 6; con- 
versely, in model GOO, where halos suffer a gradual reduc- 
tion of the available gas mass, Pop III stars continue to 
form at z £ 6, though with a declining rate. As the star for- 
mation and photoionization rate at these redshifts are well 
constrained by observations, the star formation efficiency 
and escape fraction of Pop III stars need to be lower in 
model GOO in order to match the data (see Table 1). Be- 
cause of this low /osc,iii x e»,in, reionization starts late in 
model GOO (z £ 15) and only 16% of the volume is reion- 

ized at z = 10 (reionization starts at z ~ 20 in model CF06 
and it is 85% complete by z = 10). For 6 < z < 7, QSOs, 
Pop II and Pop III give a comparable contribution to the 
total photo-ionization rate in model GOO, whereas in model 
CF06 reionization at z < 7 is driven primarily by QSOs, 
with a smaller contribution from Pop II stars only. The pre- 
dicted electron scattering optical depths are consistent with 
WMAP 3-yr data but with ~ la difference among the two 
models, leaving a chance of probing them with forthcom- 
ing CMB experiments. Of course, it would still be interest- 
ing to see whether we can find any additional detectable 
signatures, particularly with respect to the upcoming 21cm 
experiments, which is discussed in the next section. 



3 21 CM SIGNAL 

The predicted free electron fraction and gas temperature 
evolution in the redshift range 7 < z < 20 represent the 
largest difference among the two feedback models. There- 
fore, the global 21 cm signal emitted during these cosmic 
epochs might provide a viable tool to discriminate the two 
scenarios. 

The 21 cm brightness temperature, relative to the CMB, 
is given by 



5T b 



27x m (l + 5) 
(T s - r l 



0.023 



0.15 1 + zX 1 / 2 



n m h 2 io 



mK, 



(8) 



where S is the fractional overdensity, xm is the hydrogen 
neutral fraction, and T 7 is the CMB temperature at redshift 
z (see Furlanetto, Oh & Briggs 2006 for a thorough review of 
the subject). The spin temperature Ts, which represents the 
excitation temperature of the 21 cm transition, determines 
whether the signal will appear in emission (if Ts > T 7 ) or in 
absorption (if Ts < T 7 ). The evolution of the spin temper- 
ature is controlled by three competitive processes: (i) ab- 
sorption and stimulated emission of CMB photons, which 
tend to couple Ts to the CMB temperature; (ii) collisions 
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Figure 2. The best-fitting model CF06. The adopted parameters are given in Table [T] The panels show the redshift evolution of: 
(a) the volume-averaged electron and HI fraction. The arrows show an observational lower limit from QSO absorption lines at z = 6 
and upper limit from Lya emitters at z = 6.5; (b) the cosmic star formation history, with the contribution of Pop III and Pop II stars. 
Observational data are taken from the compilation of Nagamine et al. (2004 ); (c) the number of source counts above a given redshift, with 
the observational upper limit from NICMOS HUDF IIBouwens et al.ll2005l) ; (d) the electron scatte ring optical depth, with observational 
constraint s from 3-yr WM AP data; (e) Lya effective optical depth with data from lSongailal (|2004h ; (f) Lyff effective optical depth with 
data from ISongailal {2004); (g) the evolution of Lyman- limit systems with observational data from IStorrie-Lo mbardi et al.| 19941) ; (h) 
photoionization rates for neutral hydrogen, with estimates from numerical si mulations (poi nts with errorbars, [Bolton"* etaL 120051) ; (i) 
temperature evolution of the mean density IGM, with observational data from Schaye] l|l999l ). 
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Figure 3. Same as Fig. [2]but for model GOO. The adopted parameters are shown in Table IT1 
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with neutral hydrogen atoms and free electrons, which tend 
to couple Ts to the kinetic temperature of the gas Tk; (iii) 
the absorption and subsequent spontaneous emission of Lya 
photons, which mix the hyperfine levels and couples Ts to an 
effective colour temperature Tc which, in most situations, 
can be well approximated by Tk- The latter process is re- 
ferred to as Wouthuy sen- Field (W-F) mechanism and it is 
effective in the presence of substa ntial UV light, requiring 
abou t 1 Lya photon per baryon ( Wo uthuvsenl fl952 : Field 
1 19591 ). As a result, the spin temperature evolution can be 
written as 



1 s — 



1 _ J 7 



T- 1 + xcT- 1 + x a Tc 



-l 



1 + XC + Xa 



(9) 



where 



Jai R21 
T-y A21 



and 



T21 ggt 

T-y A21 



are the coupling coefficients for collisions and UV scattering, 
T21 = 0.068 K and A21 = 2.85 x 10~ 15 s _1 are the equiva- 
lent temperature and spontaneous emission rate of the 21 
cm transition, and R§i = iih k^x + n c (R21 ) is the rate 
coefficient for spin de-excitations by H-H and H-e collisions 
(UV scattering). In this an alysis, we have used the value of 
fcai tabulated in ZvecI man (20050, and the functional form 
of k2i(Tfc) reported in iLisztJ (120011) . The rate coefficient for 
the W-F mechanism is given by l|Hiratall2006l ) 

y 

where A„ = 1216 A, 7 = 50MHz is the half width at half- 
maximum of the Lya resonance, J a is the flux of Lya pho- 
tons (in cm _2 s _1 Hz _1 sr _1 ), and S a is a factor of order unity 
that accounts for spectral d istortions. We have used the nu- 
merical fits of lffiratal l|2006l ) for S a and the effective colour 
temperature Tc which appears in equation ©. 

Using equations (JSj and ([9]), we can compute the all-sky 
21 cm background signals predicted by the two radiative 
feedback models. We therefore take the redshift evolution 
of the gas kinetic temperature, Tk, and neutral hydrogen 
fraction, xhi, for model CF06 and GOO, shown respectively 
in panels (i) and (a) of Figs. [2] and [3] We compute the Lya 
background as 



J a (z) 



(1 + *) 

47T 



dz 



, dl 

dz 



I f \ — Toff (v n , Z ,z') 

-t{v, z ) e otty ' ' , 



where v = v a (l+z')/(l+z), dl/dzisthe proper line element, 
and T e fs{v a , z, z') is the effective optical depth of the IGM 
to radiation emitted at z 1 and observed at z at frequency 
v a (see section 2.2 of Salvaterra & Ferrara 2003 for a full 
description of the IGM modelling) . The comoving emissivity 
e(u, z) is computed as 



e(v,z) 



dz'L(t ZtZ ')p(z), 



2 For kinetic temperatures in the interval 1 K< Tjj <300 K we 
have used the recommended rates tabulated in column (4) of Ta- 
ble 2 and, for T > 300 K we have applied the suggested analytic 
fit. 




Figure 4. Top panel: redshift evolution of the spin (thick lines) 
and gas kinetic (thin lines) temperatures predicted by the two 
models. Solid lines refer to model CF06; dashed lines to model 
GOO. For comparison, we also show the evolution of the CMB 
temperature (dotted line). Bottom panel: corresponding evolution 
of the free electron fraction. 



where l v (t z z i) is the template specific luminosity for a stel- 
lar population of age t z z / (time elapsed between redshift z 
and z). The final value of e{v, z) is then computed by sum- 
ming over the Pop II and Pop III contribution (Salvaterra 
et al. 2006). In Fig. [4] we show the evolution of the spin tem- 
perature, gas kinetic temperature and CMB temperature for 
the two models. As expected, at very high redshifts the spin 
temperature is coupled to the CMB, because collisions are 
negligible, and it decouples from the CMB because of the 



W-F effect, which determines the evolution at 



20. As 



the star formation rate is dominated by Pop II stars, the 
Lya background is essentialy maintained by Pop II stars at 
all redshifts, and it is almost independent of the small halos 
hosting Pop III stars, which are heavily affected by radiative 
feedback. As a consequence, the spin temperature is driven 
towards the corresponding color (kinetic) temperature with 
comparable rates in the two models; the different behavior 
in the redshift range 7 z £ 19 is all due to the dissimilar 

color (kinetic) temperature evolution caused by the two ra- 
diative feedback models. In particular, in model GOO the gas 
kinetic temperature is heated above the CMB value only at 
z 15. Therefore, between 15 ^ z £ 20, Ts < T7 and we 
expect to see an absorption feature in the 21 cm signal which 
should be almost negligible in model CF06. Note that X- 
rays from supernovae, X-ray binaries and mini-QSOs can be 
an important heating agent for the neutral IGM. To cancel 
the predicted absorption feature, X-ray heating should in- 
crease the kinetic temperature from ~ 10 K to values above 
~ 30 - 40 K in the redshift (time) interval 17 < z < 20 
(~ 50 Myr). However, assuming that a fraction of 0.01 (0.1) 
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Figure 5. Predicted all-sky 21 cm brightness temperature rel- 
ative to the CMB. Solid line refers to model CF06, dashed line 
to model GOO. The upper x-axis shows a reference observational 
frequency scale. The dotted line shows an estimate of the fore- 
ground signal at these frequencies, rescaled by a factor 10 — 3 (see 
text) . 



of the total energy at Lya is emitted in X-rays (Chen & 
Miralda-Escude 2004), the temperature evolution of neutral 
regions predicted by the CF06 and GOO models is foud to 
deviate from adiabatic cooling only at z < 11 (< 13). 

Fig[S] shows the predicted 21 cm brightness temper- 
ature as a function of redshift. The upper x-axis shows 
some reference observational frequency values, computed as 
v = U2i/(1 + z), where u 2 i = 1420 MHz. The expected 
absorption signal in model GOO appears in the frequency 
range 75 - 100 MHz and has an amplitude of ~ 15 mK; 
this feature is far less evident in model CF06. At higher 
frequencies, the late reionization history predicted in model 
GOO (see the bottom panel of Fig. |4} causes a shift in the 
emission signal, which is larger than for model CF06. Note 
that since this is an all-sky signal, single-dish observations 
from existing and planned low-frequency radio-telescopes 
can reach the required mK sensitivity. The largest limita- 
tion to these measurements is the presence of strong fore- 
grounds in the relevant frequency range. The dotted line is a 
naive estimate of the contamination of a quiet region of the 
sky f rom foregrounds ( mainl y Galactic synchrotron), taken 
from iFurlanetto et all (|2006t ) . It implies that the predicted 
signals should be completely screened by foregrounds, which 
are three orders of magnitude larger. 

Since the frequency dependence of the signals and fore- 
grounds are different, the gradient of the brightness fluctua- 
tion with frequency shows some spectral features that might 
help to discriminate the signal from the relatively smooth 
foreground. Fig. |S] shows the predicted gradient for the two 
models. The largest differences are due to the deeper ab- 
sorption feature present in model GOO, and to the shift of 



Figure 6. Frequency gradient of the 21 cm brightness shown in 
Fig. \5\ adopting the same line coding. 



the emission preceding reionization. Clearly, these are of or- 
der ~ 1 — 2 mK/MHz and have to be extracted from a 
\dT B k y /dv\ ^ 3 K/MHz foreground so that their detection 
will be very challenging. A discussion of possible strategies 
for this high-precision 21 cm measurement can be found in 
IFurlanetto et all (120061). 

Following IValdes et all (|2007h . we can tentatively as- 
sume that a successful detection requires the difference in 
the brightness temperature among the two models to be 



A5T b = (5T b ) G00 - (<m) C F06 > 3 mK 
and the difference in the gradients to be, 



dST b 



dST b 



dST b 



> 0.6 mK/MHz. 



These conditions, based on the foreseen sensitivity of fu- 
ture 21 cm experiments, should take into consideration the 
difficulties associated with foregrounds removal. Using these 
limits as a guideline, Fig[7]shows that the two radiative feed- 
back models could be discriminated through the differences 
in their predicted 21 cm background signals in the observed 
frequency ranges 73-79 MHz and 82.5-97.2 MHz, which cor- 
respond to the redshift intervals 17-18.4 and 13.6-16.2, re- 
spectively. 



4 CONCLUSIONS 

Us ing a self-consistent s emi-an alytical model developed 
by IChoudhurv fc Ferraral (120051) and recently updated by 
IChoudhurv fe Ferraral (12006^ 7 we have explored the effect of 
reionization and its associated radiative feedback on galaxy 
formation. The suppression of star formation in low-mass 
galaxies due to the increase in temperature of the cosmic gas 
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Figure 7. Difference in the brightness temperature and gradients 
between the two models (see text). The shaded areas represent 
the regions where the signals will not be distinguishable given the 
foreseen sensitivity of future 21 cm experiments. Numbers along 
the curves indicate the frequency in MHz (and the corresponding 
redshift) at which the signals would be observed. 



in ionized regions has been modelled according to two differ- 
ent feedback prescriptions: (i) galaxies can form stars unim- 
peded provided that their circular velocity is larger than 
a critical threshold, which is not fixed to a constant value 
but evolves according to gas temperature (model CF06); 
(ii) depending on the mass of the galaxy, the fraction of 
gas available to star formation is reduced with respect to 
the universal value and it is fully specified by the filtering 
mass at that redshift (model GOO). We then constrain the 
four free-parameters of each model with existing observa- 
tional data, such as the redshift evolution of Lyman-limit 
absorption systems, the Gunn-Peterson and electron scat- 
tering optical depths, the cosmic star formation history, the 
number counts of high-z sources in the NICMOS HUDF, 
and thermal history of the IGM. We find that: 

(i) The two models provide equally good fit to existing ob- 
servations, which therefore are not able to break the degen- 
eracy among different radiative feedback implementations. 

(ii) Given that Pop III stars are preferentially hosted in 
the smallest galaxies at each redshift, the largest difference 
between the two models reside in the predicted Pop III star 
formation history: in models CF06, Pop III stars are fully 
expired by z ~ 6, whereas in model GOO Pop III stars con- 
tinue to form at lower z, though with decreasing rate. 

(iii) The reionization and thermal history of the IGM are 
very different: reionization starts at z £ 15 — 20, and it is 

already 85% complete by z ~ 10 in model CF06, while, at 
the same z, the ionized fraction is only 16% in model GOO. 
Both models match SDSS constraints on the evolution of 
the neutral hydrogen fraction at z < 7, but predict different 



Thomson optical depth, with r c = 0.1017, 0.0631, for model 
CF06 and GOO, respectively, in agreement with WMAP 3-yr 
data. In principle, one can use this property to distinguish 
between the two models, particularly if the observed value 
of r e is better constrained by future CMB experiments like 
PLANCK. 

(iv) Given the different gas ionization fraction and tem- 
perature evolution in the range 7 £ z ^ 20, the two models 
predict different global 21 cm background signals in the ob- 
served frequency range 75MHz u 200 MHz. The largest 

differences in the two models are represented by a ~ 15 mK 
absorption feature in the range 75-100 MHz in model GOO 
(which is nearly absent in model CF06) , and by a global shift 
of the emission feature preceding reionization towards larger 
frequencies in the same model. Single dish observations with 
existing or forthcoming low-frequency radio telescopes such 
as LOFAR, 21CMA, MWA, and SKA can achieve mK sen- 
sitivity allowing the identification of these signals provided 
that foregrounds, which are expected to be three orders of 
magnitude larger, can be accurately subtracted. 

(v) The best observational frequencies to discriminate the 
radiative feedback models through their 21 cm background 
signal are 73-79 MHz and 82.5-97.2 MHz, where the ex- 
pected differences in brightness temperatures and gradients 
are large enough to be detectable with future 21 cm experi- 
ments. 



How robust is the proposed method of differentiating 
the two radiative feedback models? It is important to stress 
that, for a given radiative feedback model, it is not possi- 
ble to find a different set of parameters which provides an 
equally good fit to the data presented in figures [2] and [3] and 
yet produce a different 21cm signal. This is because the main 
difference between the GOO and CF06 models depends on (i) 
the confinement of Pop III stars in the smallest star forming 
halos by chemical feedback, and on (ii) how radiative feed- 
back affects the gas content and therefore the star formation 
efficiency in these small mass objects. For the same reasons, 
we do not expect the resulting 21cm signals to be affected 
by varying star formation efficiencies and escape fractions 
with redshift and/or mass: at 5 < z < 20 more than 75% of 
Pop III halos are confined in a 0.5 dex mass bin centered on 
the minimum mass to form stars. Therefore we believe that 
the differences we predict in the 21cm signals are robust. 

The absorption and emission features that we find 
are comparable to the 21 cm signatures expected in mod- 
els which consider the high-redshift signals emerging from 
the Dark Age s, following the appear e nce of first lumi- 
nous sources (IChen fc Miralda-Escudi 12004 ISethil 120051 : 
Furlanctto 2006 ]), or induc ed by decaying/annihilating dark 
matter (|Valdes et al.ll2007f ). It is clear that measuring the 21 
cm background would offer valuable insights into these early 
cosmic epochs. Our analysis suggests that self-consistent 
reionization models which are compatible with a large set 
of existing observational data predict different 21 cm back- 
ground signals, which reflect how the process of star forma- 
tion in the smallest galaxies is affected by radiative feedback. 
Future 21 cm data, complemented by the available obser- 
vational data from SDSS, WMAP, thermal history of the 
IGM, high- z number counts, and the cosmic star formation 
history, might break current degeneracies and constrain the 
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reionization process and its sources over the redshift range 
6 < z < 20. 
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